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■ Abstract. A detailed description of the temporal evolution of the ther- 

modynamic fields in heavy ion collisions is presented within a hydrodynamic 
framework. Particular attention is devoted to the evolution of the collective 
^ ' flow fields and their space-time gradients. 

o 

o 

m 

O ' 1. Introduction and Motivation 
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Experimental data from RHIC , in particular the systematic analysis of anisotro- 
pics in the particle spectra, indicate that at high center of mass energies the collision 
region rapidly thermalizes and for the major part of its evolution subsequently ex- 
^ \ pands according to the laws of ideal hydrodynamics (for a recent review that collects 

the arguments and references which support this hydrodynamic scenario see 
This is an important finding by itself as apparently the generated system spends 
some time at temperatures above the anticipated QCD transition temperature and 



^ \ evolves according to the pressure, pressure-gradients and thus equation of state of 

deconfined strongly interacting matter. 

Furthermore the hydrodynamic evolution can be used to shed light even onto 
some microscopic scattering phenomena, such as signals from hard probes which 
propagate through the hydrodynamically evolving 'background' fields and radiate 
energy, as well as the late freeze-out dynamics of the reaction. Such probes are sig- 
nificantly dependent on the evolution of the macroscopic fields such as temperature 
and particle density, but also on the collective flow velocities and their gradients. 

We will here give a detailed study of the evolution of these fields. First we 
investigate a simple analytically accessible parametrization of the flow field and its 
divergence. Then we will discuss the hydrodynamic evolution of thermodynamic 
fields in a full hydrodynamic calculation to return later to the collective flow fields 
and compare to the initial analytic example. 
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2. Expansion rate and dilution rate 

Ideal hydrodynamics incorporates the continuity equations for conserved charges 
d^n^^ — where — nu^ is the current associated with the conserved charge 
whose density distribution is given by n(x) and u'^{x) is the four- velocity of the 
collective flow field. In ideal hydrodynamics the entropy density is such a conserved 
quantity. Slightly reshaping the continuity equation leads to the relation 

d^u^ = --u''d^n. (1) 
n 

The expression 9^jU^ is commonly referred to as expansion rate. Considering the 
trivial condition of vanishing collective flow = (1,0,0,0) as is the case in the 
center of the system due to symmetry, the expansion rate is simply given as = 
—h/n, the dilution rate. Assuming furthermore a power-law decay of the density n ~ 
T~" as is the case for the three-dimensional Hubble expansion of the Universe {a = 
3) or the one-dimensional longitudinal Bjorken expansion {a = 1) the expansion 
rate is given by dfj_u'^ — a/r (a is referred to as expansion parameter). 

Let us now investigate a more general situation, namely a longitudinal Bjorken 
flow field Vz — tanhyy = z/t with a transverse flow component given as jvt = 
tanh(^r). We will later see that the hydrodynamic calculations in fact suggest such 
a flow profile at freeze-out with ^ « 0.07 fm~^. It is quickly shown that this flow 
field leads to an expansion rate 

d^uf^ = T^^ cosh(^r) -t- ^ cosh(^r) -I- r^^ sinh(^r) , (2) 

which in the limit (r, z) ^ leads to 1/t + 2^ as discussed in [|^. Keeping 

the radial dependence in the expression one finds that the expansion rate increases 
with radial distance from 0.207 fm^^ at r = to 0.247 fm~^ at r = 10 fm (assuming 
T = 15 fm/c, a typical fireball lifetime). This will become important later in this 
study. 

In the case of such a non-zero transverse flow profile the spatial part of the 
four product on the right hand side of Eq. JQ) does not vanish and the expansion 
rate is no longer identically equal to the dilution rate. This becomes apparent by 
assuming a boost invariant distribution of the conserved charge (which means that 
the charge density only depends on the proper time r — \/t^ — and not on the 
space-time rapidity r\) and describing its dilution locally in terms of a power law 
such that n{r,T) ~ no{r){To/T)°' . It then quickly follows from Eq. ^ that 

rd,u^=,(a-rv,.^). (3) 
V Mr) J 

As the density is expected to drop with increasing radial distance the last term in 
the brackets is negative and gives a positive correction to a as does the preceding 
7 factor. It thus follows that in general the expansion rate (J^^u'^ is greater or equal 
to the local dilution rate a/r. 
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3. Time evolution of thermodynamic fields in central collisions 

We now investigate the evolution of the thermodynamic fields using a hydrody- 
namic model which explicitly assumes boost-invariance throughout the evolu- 
tion. The initial field configuration in the transverse plane is determined by an 
optical Glauber-parametrization and the parameters of the calculation, namely 
the equilibration time Tcqu and the constant Scqu which relates the geometrical den- 
sity to the initial entropy density are obtained by fitting to preliminary particle 
spectra from central Au-fAu collisions with ^snn = 200 GeV. Good agreement is 
achieved with Soqu = HO fm^"^ at Toqu = 0.6 fm/c [IHj, to be used in the following. 

Evolution of entropy and temperature: We first study the time evolution 
of the entropy density and the temperature at the center of the system and at 3 
and 5 fm radial distance (Fig. QJ. Initially the entropy decays with a t^^ power 
law due to the initial longitudinal Bjorken expansion (see the dashed lines). As 
time evolves, the transverse pressure gradients lead to a collective radial expansion 
and thus a transition from the initial one-dimensional to a full three dimensional 
expansion in the late stage. At late times the dashed line indicates a decay ^ t~'^. 




Fig. 1. Evolution of the entropy density (left) and the temperature (right) for 
central Au-f Au collisions at ^snn = 200 GcV. Shown are the densities and tem- 
peratures at distances i? = 0, 3 and 5 fm from the center of the system. 

The equation of state supplemented to the hydrodynamic equations of conser- 
vation of energy and momentum features a strong first order phase transition at a 
critical temperature of Tc = 164 MeV. Above that temperature the system acts like 
an ideal gas with = dp/de =1/3 and e ^ T^. For a one-dimensional Bjorken 
expansion the energy density decays as e ~ t^'^'^+^b) and thus T ~ r^^/"^. The fig- 
ure shows that this behavior is neatly followed until just before the system hits the 
mixed phase where the temperature remains unchanged while the system boils off 
its latent heat. Further cooling appears then in the resonance gas phase below Tc. 

Evolution of radial velocity: Now we consider the collective radial velocity 
which is generated through the act of the pressure gradients in the system. The 
left plot of Fig. 121 shows the radial flow velocity at distance 3 and 5 fm from the 
fireball center. As pressure gradients are larger at 5 fm than at 3 fm, the acceler- 
ation is stronger and thus the rise of Vrir). When the system reaches the mixed 
phase, the acceleration ceases and the velocity locally drops as matter continues to 
move radially outward, thus transporting the slower fluid cells from the interior to 
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larger radial distances. Only as the matter cools into the hadronic phase, pressure 
gradients pick up again and more radial flow is generated. 
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Fig. 2. Evolution of the radial velocity field. The left plot shows the velocity 

as a function of time at positions 3 and 5 fm away from the center of the system. 

The right plot shows the radial flow profile as a function of the radial distance at 4 

different times together with a linear relation w,. = with ^ = 0.07 fm~^. 

The right hand plot of Fig. [5] shows the radial flow field as a function of radial 
distance at four different times. The generated fiow field 1 fm/c after impact has an 
interesting structure, showing a non-monotonous behavior due to the mixed phase 
where no collective fiow develops. As time evolves the field rapidly adopts a linear 
behavior which is largely preserved for the rest of the evolution, and can very well 
be described by Vr{r) = ^ r tanh(^r)) with ^ = 0.07 fm~^. This 'static' fiow 
field is at the heart of the successful description of a large number of observables in 
terms of the popular 'blast wave' parametrization [{7\. 

Evolution of dilution- and expansion rates: We now return to the dilution- 
and expansion rates in the full hydrodynamic calculation. Fig. |31shows the evolution 
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Fig. 3. Evolution of the local expansion parameter a and the product of expansion 
rate and time, rd^u^^ at distances R — Q/i and 5 fm from the center. 

of the local expansion parameter a — = — ^ ^ in the left panel. It is 

clearly seen how the initial one-dimensional motion with a = 1 becomes more three- 
dimensional and just reaches a = 3 at the time of decoupling (assuming freeze-out 
at Td = 100 MeV). Note the influence of the radial position during late stages of 
the expansion. At larger radial distances the dilution rate is smaller, as supplies 
from the interior are constantly streaming outward. 
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The right panel of the figure shows the evolution of the product of expansion rate 
and time, rdf^u^. Note that with increasing the radial distance we get an opposite 
behavior as exhibited by the dilution rate. As demonstrated by the analytical 
calculation before, the expansion rate increases with radial distance. The sign of 
the difference is understood in terms of Eq. J^J. 



4. Time evolution of anisotropies in non-central collisions 

Observables from non-central collisions deliver the strongest arguments for rapid 
thermalization and a subsequent hydrodynamic expansion of the reaction region 
[ |Sj . The geometric anisotropy in non-central collision ex results in anisotropic 
pressure gradients which in turn lead to anisotropies in the collective flow field. 
Those can be characterized in terms of ep, the asymmetry of the transverse part of 
the energy momentum tensor This flow anisotropy will finally be observable as 
an azimuthal dependence of particle emission (i.e. elliptic flow) [i9|. From the left 
part of Fig. 01 which shows contours of constant energy density in the transverse 
plane at different times, it can be seen how the initial spatial eccentricity quickly 
diminishes. This becomes more quantitative in the right part of the figure where the 
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Fig. 4. Contours of the energy density in the transverse plane at four different times 
during the systems' evolution. The figure to the right shows the time evolution of 
the spatial eccentricity and the momentum anisotropy Cp, for a system with 
(RHIC 1) and without (EOS I) a phase transition. 

evolution of and Cp for a system with (RHIC 1) and without (EOS I) a first order 
phase transition is shown. The large initial geometric eccentricity leads to a rapid 
generation of momentum anisotropy, which in turn reduces the spatial eccentricity, 
and thus cuts off the further generation of momentum anisotropy. During the soft 
phase of the transition pressure gradients vanish (and thus also their anisotropies), 
which leads to an even earlier saturation of the momentum anisotropy. For systems 
at RHIC most of the anisotropy is thus generated in the plasma stage of the fireball. 
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5. Summary and Relevance 

We have given a detailed account on the evolution of the thermodynamic fields 
in systems created at RHIC within a hydrodynamic framework. This information 
is relevant for the microscopic dynamics at different levels. At early stages, hard 
probes are to be transported over the varying thermodynamic background fields to 
calculate their energy loss. At late stages, the scattering rates in the resonance gas 
needs to be compared to the macroscopic expansion and dilution rates to investi- 
gate the overall reliability of the idealized macroscopic approach which requires a 
'slow' (adiabatic) evolution and strong rescattering. Furthermore there are sensi- 
tive probes of the density around decoupling such as light nuclei, whose formation 
probability favors a large density, but their survival rate increases with decreasing 
density. Resonances reconstructed from a hadronic decay channel have similar in- 
formation content [ llOj . In summary, essentially all experimental observables are 
influenced by the temporal evolution of the thermodynamic fields presented here 
either as the signal is sensitive to all stages of the evolution and thus integrates 
the temporal history of the system or as it results from a particular stage of the 
collision with densities and their gradients giving the relevant physical background. 
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